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Abstract 

We study theoretically the renormalization of the spin-orbit coupling constant 
of two-dimensional electrons by electron-electron interactions. We demon- 
strate that, similarly to the g factor, the renormalization corresponds to the 
enhancement, although the magnitude of the enhancement is weaker than 
that for the g factor. For high electron concentrations (small interaction pa- 
rameter r s ) the enhancement factor is evaluated analytically within the static 
random phase approximation. For large r s ~ 10 we use an approximate ex- 
pression for effective electron-electron interaction, which takes into account 
the local field factor, and calculate the enhancement numerically. We also 
study the interplay between the interaction-enhanced Zeeman splitting and 

interaction-enhanced spin-orbit coupling. 
PACS numbers: 73.20.Dx, 71.45. Gm, 71.70.Ej 
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I. INTRODUCTION 



Early experimental studies of magnetotransport in two-dimensional (2D) electron 
systemsS indicated that the g factor of electrons in these systems may differ significantly 
from its bulk value. It was establishedS that the magnitude of the g factor for electrons 
confined to (100) Si surfaces exceeds g — 2 and increases from g = 2.47 to g — 3.25 with 
decreasing the concentration of electrons from 6 x 10 12 cm -2 to 10 12 cm -2 . 

Shortly after the publication of experimental resultsB it was suggested by Janaki that the 
enhancement of the g factor can be accounted for by the electron-electron interactions. The 
argument of Janak represents a 2D version of the Fermi liquid theory! and goes as follows. 
In applied weak magnetic field, B, the quasiparticle energies for the two spin projections 
can be written as 

E^(k) = E^ \k) + ^ + E^(k,E 7 (k)), 

E l (k) = ^(k)-^- + J: i (k,E l (k)), (1) 

where k is the momentum, E^°\k) = h 2 k 2 /2m is the spectrum of a free electron, Az = g^BB 
is the bare Zeeman splitting, and Y>(k,Ek) is the self-energy 

Et,i(*) = " / ?f^//(l k " k 'D/o(^ - E u (k% (2) 

where K//(?) is the Fourier component of the effective interaction between the electrons, 
and fo is the Fermi function. Solving the system Eq. ([I]), Eq. (0) in the zero-temperature 
limit, the effective g-factor can be presented as 

A* 

where 

A* z = E,(k F ) - E^kr) = x A l x . (4) 



In Eq. (f|) m* is the effective mass 



and /cp is the Fermi momentum. The enhancement factor Xz is given by 

m f 27T f 

In the random phase approximation (RPA) one hasi 

T ^ t \ 2 7re 2 ,_ N 

Veff(q) = ~ - , (7) 
e (q + V2r s k F ) 

for g < 2kp, where Sq is the dielectric constant of the material, and r s = \/2me 2 /eoh 2 kF is 

the interaction parameter of the 2D gas. With V e ff(q) in the form Eq. ([F|), m* and Xz can 
be evaluated analytically yieldingfli 

X z = J r (r s ), (8) 



™=l-—r s + r f + (l-r*)F(r s ), (9) 
m* 7r 2 

where the function Firs) is defined as 

F(r s ) = , Ts cosh" 1 (^) , r s <V2, 
nJ2 -rl V r s J 



T(r s ) = /' cos" 1 (^] , r s > y/2. (10) 



7r*/r?-2 V r * 



In the high-density limit (r s < 1) the enhancement factor (K) takes the form (see also Ref. 

Note that the theory! neglects the frequency dependence of V e ff. As a result, Eq. @ predicts 
that interactions reduce the effective mass. In fact, taking the frequency dependence into 
accountS leads to m*/m > 1 already within the RPA (see, however, the recent numerical 
simulation Jil) . 
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Later magnetotransport experimentsB0 on quantum well structures in narrow band semi- 
conductors provided an evidence for a splitting of the conduction band in a zero magnetic 
field. The analysis of the beating patterns in electron Shubnikov-de Haas oscillations led 
the authorsi'0 to the conclusion that such a splitting can be accounted for by adding the 
spin-orbit (SO) term 

H so = ak- (a xn), (12) 

to the Hamiltonian of a free electron. Here a is the SO coupling constant, k is the wave 
vector, n is the unit vector normal to the plane of the quantum well, a = (<7i, 02, 03) are 
the Pauli matrices. The term Eq. fll~2|) was first introduced by Bychkov and Rashbaffl to 
explain the experimental results on electron spin resonance!! and a cyclotron resonance of 
holeS in GaAs/AlGaAs heterostructures. 

In order to obtain more detailed information about the SO-induced splitting of the con- 
duction band, the evolution of the Shubnikov-de Haas oscillations with a tilting of magnetic 
Beld was tracedS. Subsequently, the energy spectrum of a 2D electron in a tilted magnetic 
field in the presence of the SO coupling was studied theoretically!^!!. 

Recently a zero-field splitting in different 2D systems was inferred experimentally either 
from the Shubnikov-de HaasEHI or from the commensurability oscillations^ (in a spatially 
modulated sample) patterns. 

In the domain of weak magnetic fields the Shubnikov-de Haas oscillations are smeared 
out. However, it was demonstrated both experimentally!! and theoretically!! that the SO 
coupling still manifests itself in this domain through the weak localization corrections to the 
conductance. Early worksHH in this direction were succeeded by detailed studies!!. 

In the present paper we investigate theoretically the interplay between the SO coupling 
and the electron-electron interactions. Namely, we address the question whether the inter- 
actions cause the renormalization of the coupling constant a in Eq. flT2| ) as it is the case for 
the g-factor. 
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II. RENORMALIZATION OF THE SO SPLITTING 



We assume that the bare SO coupling is weak enough ak F <C E F . With the SO term 
Eq. (|12[) the Hamiltonian of non-interacting electrons can be presented in the form 



H = E ( +\k)P + (k) + E^\k)p-(k), 
where the projection operators P + (k) and P~(k) are defined as 



(0), 



(13) 



P+(k) 



V 



p-(k) = l-P + (k), 



(14) 



so that P + (k)P (k) = 0. In Eq. (|i~4"D 0k = axctan(k y /k x ) is the azimuthal angle of the 
wave vector k. The energy spectrum consists of two branches 



Et\k) 



2m 



(15) 



Following the Fermi liquid theory, the self-energy in the presence of the SO coupling becomes 
an operator 

d 2 k' 



E(k) 



(27T) 



- 2 Veff(\k-k\) 



P+(k )f (E F - E + (k )) + p-(k )f (E F - E_(k )) 



(16) 



Our main observation is that in the presence of the interaction, the operator £(k) still 
retains the structure of Eq. ([13|) 



S(k) = S+(A:)P + (k) + 5T(fc)p-(k), 



where S ± (/c) are the scalar functions of the absolute value of the wave vector k 



(17) 



21 ' 



E± W = -2/(2^^d k - k 'D 
4/ ^2 cos (*-^')K//(|k-k'|) 

For renormalized energy spectrum, we have 



f (E F -E + (k'))+f (E F -E4k')) 
fo(E F — E + (k')) — fo(E F — E_{k')) 



(18) 



E + (k) = Ef\k) + £+(*), £L(Jfe) = P W (A;) + E~(Jfe). 



(0), 



(19) 



By solving the system Eq. ([18]) and Eq. (0), we get the following result for the renormalized 
SO splitting 

A^ = E + (k F ) - E.{k F ) = A *° , (20) 

1 — ^SO 

where Aso = lakp is the bare SO splitting and the renormalization factor is determined as 

Ago = ^h) 2 Jo ^ cos ^//(2A; F sin-). (21) 

If V e ff does not depend on (when interactions are short-ranged due e.g. to the presence of a 
gate electrode close to the 2D plane), then we have Xso — 0. However, in general, the integral 
is positive. Thus, we conclude that the exchange interaction leads to the enhancement of 
the SO coupling. Within the random phase approximation when V e ff has the form Eq. ([?]), 
the integral (|21|) can be calculated analytically and expressed through the function T{r s ) as 
follows 

^ = f- + ( 1 "^(4 (22) 

Comparison of the last expression with Eq. (|9|) indicates that 1 + Xso — rn/m*. In fact, 
this relation holds not only for V e ff (q) in the form ([7|), but for arbitrary effective interaction 
and represents a 2D version of the corresponing relation in the Fermi liquid theorjS To 
verify this relation, it is convenient to perform transformation to the real space where the 
interaction has the form V e ff(p). Then from Eqs. (|5]) and (2TD it is easy to check that 



Aso = — - 1 = T2 / dppJ?(k F p)V eff (p), (23) 
m* n Jo 



where 3\{x) is the Bessel function. Combining ( p0[) with (p3|), we get 

At 



*so 



1 + A 5 o- (24) 



A 5 o 

In Fig. 1 we plot both Xz and Xso as a function of interaction parameter r s . It is seen 
that Xz is much bigger than Xso, which has a maximum at r s = 0.52 and does not exceed 
6 percent. On the other hand, it is known that at large r s the random phase approximation 
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overestimates the screening effect which, in turn, suppresses Xso- To extend the Fermi liquid 
description to higher r s , it is customary!! to modify the random phase dielectric function! 
as follows 

e( q) =Jl- - ), (25) 

V 1 + v(q)G{q)xo J 

where v i K q)=2i\e 2 j e^q is the Fourier component of the Coulomb interaction and xo = — m /'Kh 2 
is the Lindhard susceptibility of the free electron gas. The factor G(q) (local field correction) 
describes the reduction of the screening at large q (small distances). For G(q) = we recover 
Eq. (0) for the effective interaction V e ff(q) = v(q)/e(q). 

In later works@'§ii@ a different approximation for V e ff(q) was put forward 

V eff (q) = v(q) + v\q)[l - G(q)} 2 X (q), (26) 

where x(q) is defined as 

* (?) = 1 -„(<,)(! -G W)Xa (27) 
For the local field correction G{q) the adopted the following form 

G(q) = TT=t, \ ' (28) 

lq l + q{[r s ) 



where qx{r s ) = 2a(r s )kp, and Goo(r s ), a(r s ) are the numerical factors. Eq. fl26|) is written 
neglecting the spin- fluctuation-induced vertex corrections. Combining Eqs. (|26|), (28[) and 



we get the following expression for the enhancement factor of SO coupling 

Glr s sin 0/2 



r2w ^a 2 + sin 2 0/2 

Xso = ~ — —f= I d(j) cos 



Goor, + J2{a 2 + sin 2 0/2) 



In V2 Jo (r s + sin 0/2) (a 2 + sin 2 0/2) - G^r, sin 0/2 ^/a 2 + sin 2 0/2 ' 

(29) 

In Fig. 2 we present the dependence Xso{ r s) calculated numerically within the region up to 
r s ~ 8. Followingil'IH, we took Goo(r s ) from numerical calculations (at discrete values of 
r s ) of the pair correlation function^ and following Ref. |28| assumed a(r s ) ~ l.hG^ within 



the entire domain. It is seen that instead of falling down (as in Fig. 1) Xso{ r s) increases 
with r s , when the local factor is included. Note, however, that approximately constant value 
for a(r s ) was established only within a limited interval r s < 3 in Refs. P||2"8|. We used the 
same value for calculation at higher r s in order to illustrate that the enhancement factor 
can take appreciable values in this domain. 

The alternative approach to the effective interaction in 2D electron gas with r s 3> 1 is 
described in Ref. |32|. In this paper the local field factor in the conventional form ([25|) of e(q) 
was fitted in such a way that the static characteristics of the system, calculated with this 
e(q), are consistent with the Monte Carlo results of Tanatar and Ceperley@. According to 
Ref. |32] the local field correction has the form 



Oh) = \ 



+ 



{q 2 + Ablkl) 1 / 2 (q 2 + 4b 2 2 k 2 F y/ 



(30) 



where the parameters &i(r s ), b2(r s ) are listed in Ref. [32| at discrete values of r s up to r s = 40. 
The numerical results for X so calculated for these values by substituting G(q) in the form 
(PD|) into Eq. (^) are shown in Fig. 2. They indicate that for r s ~ 10 the enhancement is 
quite pronounced. However, the reliability of both approaches can still be questioned. 



III. NON-ZERO EXTERNAL FIELD 

Now let us address the situation when the Zeeman splitting and SO coupling are present 
simultaneously. We will study the interplay between the interaction-induced enhancement 
of the g factor and of the SO coupling. First assume that Zeeman splitting is caused by a 
perpendicular magnetic field. The bare Hamiltonian in this case takes the form 



(0), 



(31) 



where the modified projection operators 



7o(fc) 



7>oi(k) = l-P + ± (k), 



(32) 



are introduced. In Eq. (|32|) 70 (A;) is denned as 

7o(*) 



T 



A z k F \ 2 A z k F 



AsnkJ A S nk 



(33) 



The bare energy spectrum is given by 

£f>(*) 



(34) 



so that the splitting of the spectrum at k = k F equals 



A = JA 2 so + A|. 



(35) 



The general expression Eq. (|16|) for the self-energy retains its form in the present case after 
changing -P ± (k) by "P^(k), where the renormalized projection operators have the form of 
Eq. © 



P+(k) 



7(*) 



l + l 2 (k) 



Pl(k) = l-P+(k), 



(36) 



I -ie^ 7(ife) y 

with renormalized parameter 7(fc), which should be determined self-consistently together 
with renormalized spectrum E±(k). Since in the present case the operators P^(k) differ 
from "P5j_(k), the consequence (|T7"D of Eq. ([IB]) is not valid anymore. Instead we get the 
following system of equations 



£ + (fc)+ 7 2 (A;)£_(A;) _ E^(k) + ^(k)E { y(k) 



l + j 2 {k) 



+ 



(37) 



d*k' K//(|k-k'|) 
(2tt) 2 l + 7 2 (fc') 



fo(E F - E + {k')) + 7 2 (fc')/o(^F - E.{k')) 



E_(k)+j 2 (k)E+(k) E ( °\k)+^{k)E^(k) 



d 2 k' K//(|k-k'|) 

(2tt) 2 1 + 7 2 (A;') 



1 + 7 2 W 



+ 



(38) 



/o(^f - £-(*')) + 7 2 (^')/o(^ - #+(*')) 



7(*) 



1 + 7 2 W 



+ 



rf 2 k 

(2r) s 



COS(0 k 



-E-(k) 
k 0K//(|k-k'|) 



7b(*0 



Ef{k)-E^{k) 



(0), 



l + 7o 2 W 

f (E F -E + (k'))-f (E F -E4k')) 



+ 



(39) 
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Subtracting Eq. ( |38| ) from Eq. (|37|), we get 



1 + 7 2 W 



E+{k)-E_{k) 



E^{k)-E^{k) 



(o). 



d 2 k 



+ i72^^' k - k ')TT^ 



f (E F -E + (k'))-f (E F -E4k')) 



+ 



(40) 



(2?r) 2 L/JVI l/ l + 7 2 (A0 

Now we can apply to Eq. and Eq. PD| ) the same argument that led to renormalization 
of Aso and Az respectively. In the zero-temperature limit this results in the following 
system of equations 



7(*) 



l+7 2 (*0 



l~7 2 (fc) 
1 + 7 2 W 



E + {k)-E_{k) 



E + {k)-E_{k) 



m 

1 X so 

m 



1 

m 



7o(fc) 
l + 7o 2 W 

l-7o 2 (fc) 
l + 7o 2 (*0 



E+\k) - E K "\k) 



Ef{k)-E^{k) 



(41) 



(42) 



Dividing Eq. ( fj2| ) by Eq. (|4T| ) we get a closed quadratic equation for 7(&) 



7 2 (fc)+ 7 (fc) 



l-7o 2 Wl-^A so 



1 = 0. 



(43) 



7o (*) 1 - 

Substituting the solution of this equation back into Eq. (p7|), we get the renormalized 
splitting of the spectrum A* = E + {k F ) — E_{k F ) 



A* _ -y (1 - 7o 2 ) 2 (l - ^A 5G ) 2 + 4 7o 2 (l - ^Xz) 2 
A 



(44) 



(l + 7o 2 )(l-^A 5 o) 2 (l-^Az) 2 * 
Using the definition (^) of 70, we can rewrite the last result in the following concise form 



A* 



A 2 



A% + A\r, VI — ^Aso 



(45) 



A \J A| + A| VI - t^Az/ ^ z t 
Finally, with the use of Eqs. (^) and (|20|) , we arrive to the conclusion that renormalized 
splitting A* is related to renormalized values A* z and AJ G in the same way as the bare 
values (Eq. (32)) 



A* = x/Af + Af. 



(46) 



Consider now the case when the Zeeman splitting is caused by a parallel magnetic field 
applied along the x-direction. Then the Hamiltonian can be written as 
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H 



h 2 k 2 

2m 



^ + iake~ i<t>k 
n 2 k 2 



2m 



(0), 



(47) 



where the energy spectrum 



±-\IA% + A| 0v ^ 



+ 2A z A so — sin0 k , 



(48) 



2m 2 

depends both on the amplitude and orientation of k with respect to the magnetic field. In 
Eq. (]47|) , the projection operators P^jj(k) are defined as 



if 1 



le 



V 



-ye 



P l(k) = l-P+(k) ! 



(49) 



/ 



with the angle ip^ related to the azimuthal angle of vector k as follows 



/ a/ccos0k 

ip k = arctan — — — s 

\ak sin 0k + — 

The bare splitting of the energy spectrum at |k| = kp is equal to 



(50) 



A(0) = J A| + A| G + 2A Z A 5G sin , 



(51) 



Performing calculations similar to those for perpendicular field, it is easy to check that in 
the present case the relation between the renormalized splitting A* and Az, Aso preserves 



the form (|5T|) 



A*(</>) = JAf + Af + 2A z A* so sin , 



(52) 



IV. CONNECTION TO THE LANDAU PARAMETERS 



The above calcualations were based on the concept of effective interaction between elec- 
trons, V e ff(q). Generally speaking, Fermi liquid theory relates the observable values to the 
bare parameters of electron gas by means of interaction functionlB having the form 



/^(k,k') = r(k,k') + (^-^')r(k,k' 



7T 



vnr 



h 2 - 



F s (k, k') + (a ■ a')F a (k,k' 



(53) 
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where a and er are spin matrices, f s (k, k') and f a (k, k') are the symmetric and anti- 
symmetric parts of the interaction function, respectively. In Eq. (|)3|) F s ( a \k, k') = 
S/^o -F* < ' a ' ) cos(/0 kk ') are dimensionless quantities. The concept of effective interaction used 
above is equivalent to the assumption f a = f s . The way to extend our theory in order to 
take into account the difference between f s and f a is to modify the self-energy Eq. ([il|) as 
follows 



' 1 w? P ^ 



V eff (\k-k\)f (E F -E + (k)) (54) 
+W eff (\k-k'\)f (E F -E_(k')) 
V eff (\k-k'\)f (E F -E4k)) 
+W eff (\k-k'\)f (E F -E + (k')) . 



Here V e ff(q) corresponds to the effective interaction in Eq. (|II]), whereas W e ff(q) accounts 
for the difference between f s and f a . It is straightforward to check that with self-energy 
operator Eq. ([54]), the projection operators P + (k) and P~{k) which, in principle, should be 
determined self-consistently, still retain the form Eq. ([T4]) . Thus we can repeat the derivation 
for the enhancement of A$o in a similar way as in Sec. II. The difference is, however, that 
the relation Eq. (p3|) does not hold anymore. Indeed, the renormalization of the effective 
mass! is now determined by the effective interaction V e ff(q) + W e ff(q) through ^- = l + ^Ff, 
while Xso is determined by the first Fourier harmonics of V e ff(q) — W e ff(q); accordingly, Xz 
is determined by the zero's Fourier component of V e ff(q) — W e ff(q). Consequently, in terms 
of the dimensionless Landau parametersi, we get the following generalization of Eq. 

At 



4 50 _ x 

1 pa • 



(55) 



A so 1 + |Ff 

Note that, while the dependence of F* on r s in two dimensions was a subject of Monte Carlo 
studiesii'i, the dependence Ff(r s ) in 2D cannot be extracted form the current literature. 
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V. CONCLUSION 



The main goal of this paper is to demonstrate that alongside with fundamental charac- 
teristics, Xz{ r s), of interacting electron gas, which describes the enhacement of the g-factor, 
and was studied in many works, there exists another fundamental characteristics Xso( r s) 
which describes the interaction-induced enhancement of the SO coupling. We calculated 
this function analytically in the limit of high concentrations and estimated numerically at 
low concentrations. Note that throughout the paper we assumed the bare SO coupling to be 
small: A$o <C Ep. However, in the limit of high concentrations (r s <C 1) the corresponding 
condition is more strict: Aso "C r s Ep, which is equivalent to a <C e 2 /eq. In the intermediate 
region r s Ep <C Aso Ep the enhancement factor is given by 

(56) 

Since experimentally the concentration of carriers is varied by changing the gate 
voltageEilBIl, there exists another simple reason for the dependence of the SO coupling 
on the concentration. Indeed, the change of the gate voltage causes the redistribution of the 
confining potential, which, in turnS, affects the parameter a. This mechanism should be 
dominant at high concentrations when Xso is small. 

Note in conclusion, that if the bare SO splitting is caused by the Dresselhaus 
mechanismil, which originates from the absence of the inversion symmetry in the bulk, the 
renormalization of the corresponding splitting of the spectrum, A^, has the same form as 
Eq. (p4|) : A* D = (1 + Xso) Ad. As a result, when both A^ and Aso are present, the splitting 
of the energy spectrum is given by the same formula as for noninteracting electronsllll 



A*(0) = yjAg + Af + 2A^A^ o sin20. (57) 

Finally, let us point out that in conventional magnetotransport oscillations experi- 
ments performed up to nowEHlZhill the typical concentrations of electrons was quite high 
~ 10 12 cm~ 2 . As a result, the typical values of the interaction parameter r s were rather 
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low (r s < 1). Only in 2D hole gasEl the condition (r s > 1) was fulfilled. For low val- 
ues of r s our theory predicts that the renormalization of the SO coupling is weak. How- 
ever, in recent experiments on the electron gases in silicon metal-oxide-semiconductor field- 
effect transistors&~0 and AlAs quantum wells0, as well as in hole gases in SiGe quan- 
tum wells@, GaAs inverted semiconductor-insulator-semiconductor structures^, and GaAs- 
AlGaAs heterostructures@ the values of r s ranged fromiH^ r s pa 6 toSB r s w 24. For 
such large r s we predict a strong renormalization of the SO coupling, which might be of 
relevance for metal-insulator transition observed in these systems. For example, the strong 
renormalization of SO coupling at large r s might cause an instability of electronic spectrum 
in a clean system, so that the system in zero magnetic field would undergo a transition 
into an exotic "chiral phase" at some critical density. To be more specific, in the latest 
publication!! the critical density for metal-insulator transition in n-type GaAs was reported 
to be n c « 1.3 X 10 10 cm -2 . To accomodate all these electrons within the lower branch of 
the spectrum, corresponding to chirality "— " (see Eq. (15)), the effective coupling constant 



a should exceed a c = I |j^r J w 3.65 x 10 12 eV-m. On the other hand, the constant a for 
relatively high-density GaAs/AlGaAs structure with n w 4.0 x 10 n cm -2 can be extracted 



from Ref. |TU| to be Aso/V^irn px 1.7 x 10~ 12 eV-m. Thus, a two times interaction-induced 
enhancement of SO coupling (although Fig.2 for r s pa 5 shows a lesser value) would make a 
chiral phase feasible. 
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FIGURES 

FIG. 1. The enhancement factors of Zeeman splitting (dotted line) and spin-orbit splitting 
(full line), calculated within the static random phase approximation, are plotted vs the interaction 
parameter r s . 

FIG. 2. The enhancement factor \so, calculated numerically with local field correction taken 
into account, is plotted vs the interaction parameter r s . Full curve corresponds to the approach 
of Refs. 4 and 29 with taken from Ref. 31 at points marked with crosses. Dashed curve is 
calculated using the local field factor taken from Ref. 32 at points marked with empty circles. 
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